Resonance enhancement of particle production during reheating 
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We found a consistent equation of reheating after inflation which shows that for small quantum 
fluctuations the frequencies of resonance are slighted different from the standard ones. Quantum 
interference is taken into account and we found that at large fluctuations the process mimics very 
well the usual parametric resonance but proceed in a different dynamical way. The analysis is made 
in a toy quantum mechanical model and we discuss further its extension to quantum field theory. 
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Inflation has become one of the most successful ideas 
of high energy cosmology today [0 . Actually most of the 
current work on experimental cosmology are based on 
ideas, or designed to test inflationary models 0. No mat- 
ter which inflationary model we assume, at the end of this 
period, it must reheat the universe. This period, called 
reheating, is very important because basically all ele- 
mentary particles populating the universe were created. 
Much of the current interest in reheating follows after the 
discovery of an exponential amplification in the number 
of the particles produced during preheating. Source of 
this exponential growth is the presence of the paramet- 
ric and stochastic resonance in the solutions of the 
equations of motion. Recent investigations address the 
problem of parametric amplification of super-horizon per- 
turbations |S|] during preheating and their consequences 
in the CMB spectrum [g. 

Because the study of models of interacting fields evolv- 
ing along the expansion of the universe is a highly non- 
trivial issue, analytical progress in this direction has been 
difficult to obtain. However, as we describe later, most 
of the main features of the observed phenomena can be 
obtained by studying a toy quantum mechanical model 
of two bi-quadratically coupled oscillators. Furthermore 
the recent interest in amplification of fluctuations can be 
described in this context. 

In this letter, we describe in a simple way how para- 
metric and stochastic resonance arise from the dynamics 
of a system of two coupled oscillators and we also iden- 
tify the main problems that occur with the standard ap- 
proaches, such as the energy conservation problem. We 
derive further in a consistent way, a semiclassical equa- 
tion of motion for the scalar condensate field R(t) (see 
Eq.(|TTl) below). From the analysis of this equation, it 
is shown that the preheating phase proceeds in a dif- 
ferent dynamical way, consistent with the Heisenberg's 
principle. This fact also leads to appropriate initial con- 
ditions - an initially non-zero value for the fluctuations 
- to allow an efficient energy transfer between the os- 
cillators. The equations obtained (Eq.([To|) and Eq.jTl])) 
can be considered as a classical system, allowing us to 
perform a complete numerical study |1C| . A first numer- 



ical study shows a novel saturation effect in the small 
i?-region of the equation for the fluctuation. This effect 
holds at the beginning of the energy transfer between the 
oscillators. We extend our results to quantum field the- 
ory, and take into account the expansion of the universe. 
This is straightforward, and the results obtained in our 
previous analysis can be extended with slight modifica- 
tions. 

The model we consider is a classical oscillator X in- 
teracting with a purely quantum mechanical oscillator Y 
described by the Lagrangian 
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g 2 X 2 Y 2 , (1) 



where X{t) corresponds to the classical inflaton field <fi(t) 
and Y(t) the created scalar field x(x) of inflationary cos- 
mology respectively. For reasons of simplicity we will 
first neglect the expansion of the universe, and discuss at 
the end the effects and possible modifications of includ- 
ing such a friction term. The equations of motion are 
then corresponding to two coupled harmonic oscillators 
with time-dependent frequencies 



X + {u} 2 + g 2 Y 2 )X = 0, 



Y +{Q 2 +g 2 X 2 )Y = 



(2) 



(3) 



Let us now introduce time-independent creation and an- 
nihilation Heisenberg operators a' and a respectively us- 
ing the Ansatz 



Y = f{t)a + t{t)a\ 



(4) 



where, due to the standard commutator [a, a^] = 1 for 
creation-annihilation operators, f(t) satisfies the Wron- 
skian condition 

/*(*)/(*) - /*(*)/(*) = (5) 

Inserting the Ansatz of Eq. (||) into Eq. (|J) , we find that 
f(t) satisfies the equation of motion 
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f(t) + (tf+g 2 X 2 (t))f(t)=Q, 



(6) 



where the normalization is fixed by the condition (||). 
If we consider Y in the frame of quantum field theory, 
Eq. (|^) would be the equation for the zero-momentum 
term of its Fourier expansion. Motivated by the plane 
wave expansion of a quantized-field operator we use the 
following Ansatz for the function-coefficient f(t) 
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exp(-i / W(t')dt' 



(J) 



which satisfies automatically Eq. (|5|). Inserting this 
Ansatz into Eq. (0), we obtain the nonlinear differen- 
tial equation for W(t) 




+ W 2 = m 2 (i), 



(8) 



where m 2 (t) = Q 2 + g 2 X 2 (t). Eq. (|8|) describes a 
non-linear dissipative differential equation for W{t). In- 
stead of solving this rather cumbersome equation nu- 
merically, we will consider the Hamiltonian associated 
to the Lagrangian of Eq. ([!]) and construct a classi- 
cal effective Hamiltonian by projecting onto the sub- 
space of the number operator h = a^a, generated by 
{|4) , 1 3) , 1 2) , |1) , |0)}, which we will call S 5 . The vac- 
uum is as usual defined by the condition a |0) =0. The 
choice of the non-trivial subspace S5 is two fold. First, 
it attempts a fully quantum-mechanical treatment of this 
problem, allowing the possibility of quantum interference 
and other quantum effects in the effective classical the- 
ory, and secondly it allows an exact analitical solution 
of the quantum sector of the model. We first define the 



classical variable R 2 (t) = 
reads 



Y< 



0). The operator Y 2 
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and can be diagonalized in S5 explicitly Y 2 Vi — 

L J nrn 

AjVi, where Vi and \i are the eigenvectors and eigenvalues 
respectively. Because of the similar structure of both 

; 2 

operators, Y 2 and P Y — Y , we can use the Ansatz of Eq. 
([|), and the effective Hamiltonian can be diagonalized 
exactly as 
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[ p x +Pr + ^2+ " * + m 2 {t)R 2 ) yi 5, 



(9) 



where i, j — 0, ..4. Up to irrelevant constants, this Hamil- 
tonian describes the effective dynamics of two classical 



variables X and R, which corresponds to the original 
classical and quantum oscillators. From Eq. the 
equations of motion are 



X + (lu 2 + g 2 R 2 )X = 0, 



R+{n 2 +g 2 X 2 {t))R~ T ^ = 



(10) 



(11) 



The above equations actually correspond to two classical 
coupled particles-like degrees of freedom. Note that the 
centrifugal term keeps the fluctuation R away from zero, 
consistent with Heisenberg's principle and giving us ap- 
propriated natural initial conditions for numerical studies 
p0| . These equations are very simple and represent the 
starting point of our analysis. 

A similar derivation of the classical Hamilton dynamics 
of a classical oscillator interacting with a quantum me- 
chanical oscillator has been made in jflj , in the context of 
scmiquantum chaos. Nevertheless, the chosen classical- 
variable parametrization leads to cumbersome nonlinear 
dissipative differential equations. 

We note here, that backreaction is taken into account 
in Eq. (|To|), through the dependence on the classical 
condensate variable R. We start the discussion of the 
classical dynamics by considering the perturbative regime 
of the equations of motion, at small coupling constant 
g 2 < 1. For g 2 = 0, a solution of Eq. @ is XJt) = 
Xos'm(ojt). Inserting this expression into Eq. (J6|) we 
obtain a Mathieu equation 



f"(z) + [A(q)-2qco S 2z]f{z)=Q, 



(12) 



where z = cut, q = g 2 Xl/Auj 2 and A(q) = fl 2 /uj 2 + 2q. 
This equation is well-known and has been used to dis- 
cuss different mechanisms of particle production during 
the reheating period of the universe, such as paramet- 
ric and stochastic resonance ||- ||. In fact, for 2q < A 
this equation leads to narrow parametric resonance pj j 
and for 2q > A leads to broad parametric resonance |4 ] 
(called stochastic resonance if expansion of the universe 
is considered). From the definitions above for A and q, 
we observe that only the first region is allowed. The limit 
curve A = 2q in the space of parameters can be discussed 
by considering the limit case fi 2 = of the equivalent 
equation 



f{t) + (n 2 + g 2 X 2 sin 2 (ujt))f(t) = 0. 



(13) 



This equation leads also to a linear growth of the ampli- 
tude f(t), close to the minimum of the potential for X{t). 
This can be explained intuitively because this equation 
reads f(t) = at these points (see Fig. 1), and hence 
has solutions which raise linearly with time. 

However, it is natural to ask whether or not a non- 
exactly oscillatory form of X will lead to an exponential 
growth of f(t). Moreover, we know from this simple sys- 
tem that as soon as f(t) start to growth, the functional 
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FIG. 1. We show the oscillations of the source frequency 
(see Eq. (|l3|)) which is responsible for parametric resonance 
(curve 1) and the response f(t) (curve 2). Note that when the 
X oscillations pass through zero, the maximun amplifications 
occur. 
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FIG. 2. We plot the potential of Eq. (|l^) in a convenient 
scale. Note the the potential has a wall close to origen R = 0. 




FIG. 3. This plot shows the parametric amplification of 
R. In the context of a quantum field theory, this result is 
interpreted as a particle production process. 




FIG. 4. In the small amplitudes region, close to the res- 
onance conditions, there is a novel saturation effect due to 
detuning coming from nonlinearities of the equation of mo- 
tion. See explanation below Eq. dli|). 



form of X will change, not only in amplitude but also in 
frequency. Also, attention must be down to solutions of 
Eq. ( fTl| ) and not of Eq. (^j), because in general its solu- 
tions does not satisfy the Wronskian condition Eq. (||). 
The best way to take into account the en erg y conserva- 
tion, is by studying a numerical solution jl(J . Even, the 
system ([lO jll ) can be seen as a classical system allowing 
us to perform in a easier way such a study. 

So, let us return to equations ( |To| ) and (pi]). The sec- 
ond equation can be interpreted as a particle moving in 
the presence of the potential 



V(R) = \m\t)I? + ±- 2 . 



(14) 



Figure 2 shows a plot of the potential for a convenient 
value of m 2 . We see that R = acts as an infinite wall, 
which divides the coordinate space into two disconnected 
regions. Depending on which sign of the root we choose, 
the particle will remain in one of these two regions. 



In order to analyze this potential, there are two scales 
to be considered. For large amplitudes R l/\/2m, the 
potential looks essentially like a parabola with an infi- 
nite wall at R = (see Fig. 2). The solution of the 
equation without the non-linear term, which leads to the 
infinite wall, has resonances at fl n = 2lo /rt, as in the 
conventional parametric resonance (see [|12[). One could 
naively expect to find resonances of the complete equa- 
tion at twice this value because of the presence of the wall 
at R = 0. But this is not true. The reason for this be- 
havior is as follows. The solution of the problem without 
the wall-term is invariant under the shift ut — > Lut + n. 
But this corresponds to the same solution of the system 
in the presence of the wall, when the phase is shifted by 
7r as R approaches zero. Figure 3 shows the growth of 
the amplitude of oscillations for the main resonance band 
il = 2lu. This curve has a similar shape and properties of 
the one corresponding to narrow parametric resonance, 
Fig. 2 in §. 
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x" + 2H X ' + a 2 (M 2 +g 2 ( f> 2 )x = : 



(17) 




FIG. 5. This plot shows how R ceases to grow exponen- 
tially after certain time, when the expansion of the universe 
is considered. The resonance structure of the secondary max- 
ima is given by the detuning effect due to the presence of the 
a 2 -factor in the effective frequency term of Eq. (fell) . 



For small amplitudes R ~ l/y/2m, we can expand the 
potential around the time dependent minimum R m {t) = 
±l/y/2m(t) 



V(R) ~ -m 2 it) + 2{R ~ R m ) 2 m 2 (t). 



(15) 



From (|l5|), we see that for amplitude values close to the 
minimum, the potential develops a linear term, which 
leads to non-linear oscillations. Such oscillations hold 
in an almost energy conserved environment, because the 
amplitude of R is small, so this effect is a genuine one 
and holds during the beginning of the transfer of energy. 
Fig. 4 shows the typical behavior of these oscillations. 
There are two regions separated by the size of the am- 
plitude. When the amplitude is small and a resonance 
condition is fulfill, the linear term drives the system to 
larger values of the amplitude, which is characteristic of 
the broad parametric resonance (see Fig. 3 in ||). Nev- 
ertheless, when the amplitude grows the non-linear term 
of the equation becomes important and breaks the reso- 
nance tuning and the growth of the amplitude saturates, 
while the amplitude starts to decrease. Most of the re- 
sults obtained above can be still applied when we consider 
the expansion of the universe. 

In quantum field theory the process follows similar 
lines. If we consider the inflaton field cf> and the created 
field x coupled minimally to gravity, the Lagrangian can 
be written as a sum C = Co(<f>) + £q(x) + £i{4>, X)i where 



2 J.21 



(16) 



where the interaction is Cj — — ^g 2 (/> 2 X 2 - Let us con- 
sider <fi(t) and x{t) as homogeneous fields, but the former 
as a classical and the later as a quantum operator field. 
Assuming a flat FRW metric, the equations of motion 
obtained from the Lagrangian (fLq) are 



b" + 2H^' + a 2 (m 2 + g\ 2 )4> = Q, 



(18) 



where a is the scale factor, a prime ' means a derivative 
respect to conformal time rj = J dt/a(t), and TL — a' /a. 
Defining the conformal fields as X = a<j> and Y = ax, we 
find the equations 



X" + {a 2 m%+g 2 Y 2 - — )X = 0, 
v a 



Y" + (a 2 M 2 + g 2 X 2 )Y = 0. 

x a 



(19) 



(20) 



If we neglect the coupling between x with <f> and during 
the oscillations of the inflaton, the universe expands as 
a one dominated by matter. Under such circunstances 
the last term inside the brackets can be dropped out. 
Moreover, we have checked numerically that the last term 
a" J a can be neglected, even if the coupling is considered. 
The semiclassical approximation leads to the equation 



R" + (a 2 M 2 + g 2 X 2 



)R 



1 

AR? 



0. 



(21) 



and similarly the corresponding for X (Eq.(19)) but re- 
placing Y 2 by R 2 . If we consider the corrections beyond 
the long-wavelenght limit of the field Y, an infinite set 
of non-linear coupled differential equations for the modes 
Rk prevents an analogous explicit description as the one 
given by Eq. (pl|). Nevertheless, this observation allows 
a perturbative treatment of this problem by a gradient 
expansion of the fields, but this is beyond the scope of 
our analysis. Because the fields are scaled by a, the ex- 
pansion of the universe should damp the oscillations of 
the created field R (or x) in Fig- 1 and Fig. 3. In Fig. 5 
the combined effects of backreaction and the expansion 
of the universe in the evolution of the field R is shown. 
As a consequence, R ceases to grow exponentially, after 
certain time. The resonance structure of the secondary 
maxima can be explained by the detuning effect due to 
the presence of the scale factor in the effective frequency 
term of Eq. @. 

In this letter we have described a model to study re- 
heating after inflation. We have showed that a simple 
model of two coupled oscillators can reproduce the main 
resonances effects, responsible for particle creation during 
preheating. In addition we have derived a semiclassical 
equation describing particle production consistent with 
Heisenberg's uncertainty principle. From this equation 
we have also found that at small amplitude of the oscil- 
lations, i. e. when the energy transfer starts, there is 
no resonance amplification. This seems to indicate that 
there is an energy threshold to get efficient particle pro- 
duction. 
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